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Abstract We present recent Monte Carlo results on surfaces of bcc-structured binary 
alloys which undergo an order-disorder phase transformation in the bulk. In 
particular, we discuss surface order and surface induced disorder at the bulk 
transition between the ordered (DO3) phase and the disordered (A2) phase. 
An intricate interplay between different ordering and segregation phenomena 
leads to a complex surface behavior, which depends on the orientation of the 
surface under consideration. 

The structure and composition of alloys at external surfaces and internal 
interfaces often differs significantly from that in the bulk. In most cases, this 
refers only to very few top layers at the surface, over a thickness of order 1 
nm. In the vicinity of a bulk phase transition, however, the thickness of the 
altered surface region can grow to reach mesoscopic dimensions, of order 10- 
100 nm. If the bulk transition is second order, for example, the thickness of 
the surface region is controlled by the bulk correlation length, which diverges 
close to the critical point [1]. Close to first order bulk transitions, mesoscopic 
wetting layers may form [2]. 

While these various surface phenomena are fairly well understood in simple 
systems, such as surfaces of liquid mixtures against the wall of a container, 
the situation in alloys is complicated due to the interplay between the local 
structure, the order and the composition profiles. In alloys which undergo 
an order /disorder transition, for example, the surface segregation of one alloy 
component can induce surface order[3, 4, 5] or partial surface order[6, 7] at 
surfaces which are less symmetric than the bulk lattice with respect to the 
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ordered phase. Even more subtle effects can lead to surface order at fully 
symmetric surfaces[8, 9, 10, 11]. Furthermore, different types of order may be 
present in such alloys[12], which can interact in a way to affect the wetting 
behavior significantly [14, 15]. 

In this contribution, we discuss a situation where such an interplay of seg- 
regation and different types of ordering leads to a rather intriguing surface 
behavior: Surface induced disorder in a binary (AB) alloy on a body centered 
cubic (bcc) lattice close to the first order bulk transition between the ordered 
DO3 phase and the disordered bulk phase. Surface induced disorder is a wet- 
ting phenomenon, which can be observed when the bulk is ordered and the 
surface reduces the degree of ordering - usually due to the reduced number 
of interacting neighbors[13, 16, 17]. A disordered layer may then nucleate at 
the surface, which grows logarithmically as the bulk transition is approached. 
According to the theoretical picture, the surface behavior is driven by the de- 
pinning of the interface between the disordered surface layer and the ordered 
bulk. 

In order to study the validity of this picture, we consider a very idealized 
minimal model of a bcc alloy with a DO3 phase: The alloy is mapped on an 
Ising model on the bcc lattice with negative nearest and next nearest neighbor 
interactions. The Hamiltonian of the system then reads 

7i = Vj2 SiSj + aVj2 S i S i -HJ2 S *> (!) 

(ij) ((ij)) » 

where Ising variables S = 1 represent A atoms, S = — 1 B atoms, the sum (ij) 
runs over nearest neighbor pairs, ((ij)} over next nearest neighbor pairs, and 
the field H is the appropriate combination of chemical potentials fj,A and \xb 
driving the total concentration c of A in the alloy, (c = ((S) + l)/2). 

The phases exhibited by this model are shown in Figure 1 In the disordered 
(A2) phase, the A and B particles are distributed evenly among all lattice 
sites. In the ordered B2 and DO3 phases, they arrange themselves as to form 
a superlattice on the bcc lattice. The parameter a was chosen a = 0.457, such 
that the highest temperature at which a DO3 phase can still exist is roughly 
half the highest temperature of the B2 phase, like in the experimental case of 
FeAl. The resulting phase diagram is shown in Figure 2. 

In order to characterize the ordered phases, it is useful to divide the bcc 
lattice into four face centered cubic (fee) sublattices as indicated on Figure 1, 
and to define the order parameters 

^1 = ((S) a + (S) b -(S) c -(S) d )/2 

4>2 = ((S)a-{S) b + (S) c -(S) d )/2 (2) 

V>3 = ((S) a -(S) b -(S) c + (S) d )/2, 

where {S) a is the average spin on the sublattice a. In the disordered phase, 
all sublattice compositions are equal and all order parameters vanish as a 
consequence. The B2 phase is characterized by tpi / and the DO3 phase by 
V>i / and tp2 = ±V"3 / 0- The two dimensional vector (^2,^3) is thus an 
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Figure 1. H/V 

Figure 2. 

Phase diagram in the T—H plane. First or- 
Ordered phases on the bcc lattice: (a) dis- der transitions are solid lines, second order 
ordered A2 structure, (b) ordered B2 and transitions dashed lines. Arrows indicate 
(c) DO3 structure. Also shown is assign- positions of a critical end point (cep) and 
ment of sublattices a, b, c and d. a tricritical point (tcp). 



order parameter for DO3 ordering, and the latter can be characterized conve- 
niently in terms of its absolute value 

V>23 = \j(lp2 2 + ^3 2 )A (3) 

We have studied free (110) and (100) surfaces of this model at the temper- 
ature T = 1 ksT/V in the DO3 phase close to the transition to the disordered 
phase. To this end, we have first located the transition point very accurately 
by thermodynamic integration [18], Hq/V = 10.00771 [1]. We have then per- 
formed extensive Monte Carlo simulations of slabs each 100-200 layers thick, 
with free boundary conditions at the two confining (110) or (100) planes, and 
periodic boundary conditions in the remaining directions. 

In all of our simulations, the average value of the Ising variable in the 
top layer was one, i.e., the top layer was completely filled with A atoms. 
Having stated this, we shall disregard this layer in the following and discuss 
the structure starting from the next layer underneath the surface. The layer 
order parameters ipi(n) and the layer compositions c(n) can be determined 
in a straightforward manner for (110) layers, since they contain sites from all 
sublattices. In the case of the (100) layers, it is useful to define c and the ipi 
based on the sublattice occupancies on two subsequent layers. 



Fig. 3 shows the calculated profiles for two choices of H close to the tran- 
sition. One clearly observes the formation of a disordered film at the surface, 
which increases in thickness as the transition point is approached. The film is 
characterized by low order parameters ipi and ^23 , and by a slightly increased 
concentration c of A sites. The structure very close to the surface depends 
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on its orientation: The composition profiles display some characteristic oscil- 
lations at a (110) surface, and grow monotonously at a (100) surface. The 
order ip23 drops to zero. The order ipi drops to zero at the (110) surface, and 
at the (100) surface, it changes sign and increases again in the outmost two 
layers. The latter is precisely an example of the segregation induced ordering 
mentioned earlier [3, 4, 5]. 



(110) surface (100) surface 




Layer Layer 



Figure 3. 

Profiles of the total concentration (top) and order parameters tpi (bottom, circles) and t[>23 
(bottom, squares) at the (110) surface (left) and at the (100) surface (right) for different 
fields H in units of V as indicated. 



We can thus distinguish between two interesting regions in these profiles: 
The near-surface region, where the properties of the profiles still reflect the 
peculiarities of the surface, and the interfacial region, where the profiles are de- 
termined from the properties of the interface separating the disordered surface 
film from the ordered bulk. 

The structure of the profiles in the interfacial region is basically determined 
by the fluctuations of the interface, which are characterized by a transverse 
correlation length £m . The latter is in turn driven by the thickness of the film 
and the interfacial tension a or, more precisely, by a rescaled dimensionless 
interfacial tension 

1/w = 4^ b 2 a/k B T, (4) 

with the bulk correlation length £5. The renormalization group theory of 
critical wetting[2], which should apply here[16], predicts that the transverse 
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correlation length diverges according to a power law 

i\\<x^={H Q -H)- v W (5) 

as Hq is approached, with the exponent i/y = 1/2. With this knowledge, 
one can calculate the effective width of the order parameter profiles, £j_ oc 
sj — Luln(Ho — H), the profiles of layer susceptibilities etc. We have examined 
these carefully at the (110) interface and the (100) interface, both for the order 
parameters ipi and ip23 , and we could fit everything nicely into the theoretical 
picture. Our further discussion here shall focus on the near-surface region. 

Assuming that the order in the near-surface region is still determined by 
the fluctuations of the interface, the theory of critical wetting predicts a power 
law behavior 

lPa ;1 OC (Ho - H)^ (6) 

for value if))a, 1 of the order parameter ip a directly at the surface, regardless 
of the structure of the surface. Figure 4 shows that, indeed, the surface order 
parameter ^23,1 decays according to a power law at both the (110) and the 
(100) surface, with the same exponent and the exponent is in both cases 
identical within the error, [3\ = 0.618. Furthermore, we notice strong finite 
size effects close to the Hq. Since these are asymptotically driven by the ratio 
(L/£||), they can be exploited to determine the behavior of £y as the phase 
transition is approached, i.e., the exponent i/y. The finite size scaling analysis 
yields u» = 1/2 [19], in agreement with the theory. 



(110) surface 



(100) surface 















Slope: 0.618[4] 
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Figure 4. 

Order parameter ^23 at the surface at the (110) surface (left) and the (100) surface (right) 
vs. (Ho — H)/V for different system sizes L x L x D as indicated. Solid line shows power 
law with exponent /3i = 0.618. 



The profiles of -023 thus seem entirely determined by the depinning of the 
interface, in agreement with the standard theory of critical wetting. The 
situation is however different when one looks at the other order parameter, 
ipi. This is not particularly surprising in the case of the (100) surface. We 
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have already noted that this surface breaks the symmetry with respect to 
the ipi ordering, hence the segregation of A particles to the top layer induces 
additional -01 order at the surface (Figure 5). The (110) surface, on the other 
hand, is not symmetry breaking. The order ipi decays at the surface, yet with 
an exponent (3\ = 0.801 which differs from that observed for ^23 (Figure 6). 
Even more unexpected, the finite size effects cannot be analyzed consistently 
with the assumption that the transverse correlation length diverges with the 
exponent i/y = 1/2, but rather suggest i/y = 0.7 ± 0.05. The order parameter 
fluctuations of tp\ at the surface seem to be driven by a length scale which 
diverges at Hq with an exponent different from that given by the capillary 
wave fluctuations of the depinning interface. 




0.001 0.01 

(H„-H)/V 

Figure 5. 

Order parameter tpi in the first layers un- 
derneath the (100) surface vs. (Hq — H)/V 
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Figure 6. 

Order parameter ip\ at the (110) surface 
vs. (Ho — H)/V for different system sizes 
L x L x D as indicated. Solid line shows 
power law with exponent f3\ = 0.801. 



To summarize, we have seen that the phenomenology of surface induced 
disorder in a relatively simple bcc alloy with just two coupled types of order- 
ing is much more complex than predicted by the standard theory of surface 
induced disorder and critical wetting. Looking at our profiles, we were able 
to distinguish between two regions, the near-surface region and the interfacial 
regions. In the situations studied in our simulations, it seemed that these 
regions could be well separated from each other. Some rapid variations of 
the profiles in the near-surface region, are followed by smooth changes in the 
interfacial region. The local surface structure affects the total composition 
profile relatively strongly, and has practically no influence on the profile of 
the DO3 order, ^23 • In the case of the B2 order, ipi, the situation is more 
intriguing. The symmetry breaking (100) surface induces local order in the 
near-surface region which apparently does not couple to the interface. At 
the non-symmetry breaking (110) surface, ipi was found to exhibit qualita- 
tively new and unexpected power law behavior as the wetting transition is 
approached. 

The last observation clearly requires further exploration in the future. The 
picture will be even more complex in situations where the near-surface profiles 
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and the interfacial profiles cannot be separated any more. We expect that this 
could be the case, e.g., at (111) surfaces, which break the symmetry with 
respect to both B2 and DO3 ordering. 
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References 

[1] For reviews on surface critical phenomena see K. Binder in Phase Tran- 
sitions and Critical Phenomena, Vol. 8, p. 1 (1983), C. Domb and J.L 
Lebowitz eds., Academic Press, London; S. Diehl, ibid, Vol. 10, p.l. (1986) 

[2] For reviews on wetting see, e.g., P. G. de Gennes, Rev. Mod. Phys. 57 
827 (1985); S. Dietrich in Phase Transitions and Critical Phenomena, C 
Domb and J.L. Lebowitz eds (Academic Press, New York, 1988), Vol 
12; M. Schick in Les Houches, Session XLVIII - Liquids at Interfaces, J 
Charvolin, J. F. Joanny, and J. Zinn-Justin eds (Elsevier Science Publish- 
ers B.V., 1990). 

3] F. Schmid, Zeitschr. f. Phys. B 91, 77 (1993V 

4 S. Krimmel, W. Donner, B. Nickel, and H. Dosch, Phys. Rev. Lett. 78, 
3880 (1997). 

[5] A. Drewitz, R. Leidl, T. W. Burkhardt, and H. W. Diehl, Phys. Rev. Lett. 
78, 1090 (1997); R. Leidl and H. W. Diehl, Phys. Rev. B 57, 1908 (1998); 
R. Leidl, A. Drewitz, and H. W. Diehl, Int. Journal of Thermophysics 19, 
1219 (1998). 

[6] H. Reichert," P.J. Eng, H. Dosch, I.K. Robinson, Phys. Rev. Lett. 74, 2006 
(1995). 

[7] F. Schmid, in Stability of Materials, p 173, A. Gonis et al eds., (Plenum 

Press, New York, 1996). 
[8] L. Mailander, H. Dosch, J. Peisl, R.L. Johnson, Phys. Rev. Lett. 64, 2527 

(1990). 

[9] W. Schweika, K. Binder, and D. P. Landau, Phys. Rev. Lett. 65, 3321 
(1990). 

10] W. Schweika. D.P Landau, and K. Binder, Phys. Rev. B 53, 8937 (1996). 
11] W. Schweika, D.P. Landau, in Computer Simulation Studies in 

Condensed-Matter Physics X, P. 186 (1997). 
[12] for reviews see D. De Fontaine in Solid State Physics 34, p 73 H. Ehren- 
reich, F. Seitz and D. Turnbell eds., Academic Press, New York 1979; K. 
Binder in Festkorperprobleme (Advances in Solid State Physics) 26, p 133, 
P. Grosse ed., Vieweg, Braunschweig 1986. 
R. Lipowsky, J. Appl. Phys. 55, 2485 (1984). 

E. H. Hauge, Phys. Rev. B 33, 3323 (1985). 

D. M. Kroll, G. Gompper, Phys. Rev. B 36, 7078 (1987); G. Gompper, 
D. M. Kroll, Phys. Rev. B 38, 459 (1988). 
D. M. Kroll and R. Lipowsky, Phys. Rev. B 28, 6435 (1983] . 
H. Dosch, Critical Phenomena at Surfaces and Interfaces (Evanescent X- 
ray and Neutron Scattering), Springer Tracts in Modern Physics Vol 126 
(Springer, Berlin, 1992). 
K. Binder, Z. Phys. B 45, 61 (1981). 

F. F. Haas, Dissertation Universitat Mainz (1998); F. F. Haas, F. Schmid, 
and K. Binder, in preparation. 



[18 
19 



